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FOREWORD 


The work reported herein was conducted by the Environmental 
Research Institute of Michigan for Supporting Research, one of eight 
projects of AgRISTARS, the program for Agriculture and Resources 
Inventory Surveys through Aerospace Remote Sensing. AgRISTARS is a 
six-year program of research, development, evaluation, and application 
of aerospace remote sensing for agricultural resources which was 
initiated in Fiscal Year 1980. AgRISTARS is a cooperative effort of 
five federal agencies of the United States — Department of Agriculture, 
USDA; National Aeronautics and Space Administration, NASA; Department 
of Commerce, USDC; Department of Interior, USDI; and the Agency for 
International Development, USAID. 

The goal of the program is to determine the usefulness, cost, and 
extent to which aerospace remote sensing data can be integrated into 
existing or future USDA systems to improve the objectivity, reliability, 
timeliness, and adequacy of information required to carry out USDA 
missions. The overall approach Is comprised of a balanced program of 
remote sensing research, development, and testing which addresses 
domestic resource management as well as commodity production infor- 
mation needs in both domestic and foreign regions. 

The Supporting Research Project's objective is to develop improved 
technology through research to aseess crop area and condition, in 
support of other projects which address improved commodity production 
forecasts in foreign and domestic areas, early warning of changes in 
production conditions and expectations, improved yield estimation, and 
other applications. Dr. Jon D. Erickson, NASA Johnson Space Center, 

Code SF3, is the NASA manager of the Supporting Research Project and 
Thomas Pendleton, Code SF3, the Technical Coordinator for the reported 
effort. 
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This work was carried out within the Infrared and Optics Division, 
headed by Richard R, Legault, a Vice-President of ERIK, under the 
technical direction of Robert Horvath, Program Manager, and William A, 
Malila, Leader of the Objective Labeling Technology Development Task, 
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INTRODUCTION 


The cyclic coverage of satellite-borne remote sensing systems, 
such as Landsat, and the availability of spatially registered data have 
allowed the utilization of temporal information in analysis and dis- 
crimination of crop typeB. While the temporal-spectral pattern of 
development has proven to be an important piece of information, par- 
ticularly in analyst identification of the crop type of data samples 
[1], its use has been somewhat hindered by the intermittent nature of 
the Landsat observations and confusion caused by differences in the 
planting dates and development patterns of neighboring fields of the 
same crop. This shortcoming may be overcome by estimating the overall 
pattern of spectral development for each field or pixel from the set 
of samples provided by Landsat. Based on this more complete description 
of sample characteristics, important information can be extracted for 
use in analysis and classification. This report presents an algorithm 
designed to estimate one important piece of information — the relative 
crop calendar shift — from observations of spring small grains [2,3]. 
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BASIS OF CROP CALENDAR SHIFT ESTIMATION 
FROM LANDSAT SPECTRAL DATA 

Temporal-spectral profile characterization plays a key role in the 
utilization of multitemporal Landsat data for crop calendar shift 
estimation. This section reviews the agrophysical basis of such 
characterization, describes the crop calendar shift estimation concept, 
and discusses profile modeling of spring small grains. 

2.1 TEMPORAL-SPECTRAL PROFILE CHARACTERIZATION 

The sequence of development stages through which plants of a given 
crop type progress in the course of a growing season is accompanied by 
a related sequence of spectral reflectance stages. While not every 
stage of physiological development will necessarily result in a dif- 
ferentiable set of reflectance properties, those stages corresponding 
to observable differences in plant morphology or canopy structure will 
influence reflectance. Since most of the morphological changes occur 
gradually rather chan abruptly, particularly when viewed at the popu- 
lation (field or pixel, i.e., picture element) level, it is reasonable 
to assume that the spectral development of a field could be represented 
by a continuous function. 

The intermittent observations from a remote sensing platform such 
as Landsat can be viewed as discrete samples from such continuous 
functions, and can provide important information for crop identification 
and condition assessment. However, if the continuous function repre- 
sented by the Landsat samples could be reconstructed, considerably more 
complete information could be available. Profile characterization is 
the process of fitting mathematical functions (profiles) to a set of 
observations of a given scene element. 
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While the temporal pattern of any spectral band or variable could 
be fit with a mathematical function, those spectral features most closely 
correlated to describable plant or field phenomena are probably of 
greatest use. In particular, indicators of "greenness" (spectral 
variables which are sensitive to the green vegetation component of 
agricultural signatures) such as the Greenness variable produced by the 
Tasseled-Cap transformation [4,5] often exhibit temporal patterns which 
are relatively easy to characterize mathematically and contain sub- 
tantial information. 

Most crops follow a general pattern consisting of a period of in- 
creasing greenness to a maximum value, followed by a period of declining 
greenness. The rate of greenness increase or decline is closely tied 
to the rates of change in leaf area and canopy closure, and the rates 
of plant maturation. Specifically, we would expect a low rate of 
greenness increase as the plants emerge and begin to affect the spectral 
reflectance properties of the field, followed by a more rapid rate of 
increase corresponding to the period of most rapid bio-mass production. 
Then, as the maximum greenness is approached, increases in greenness 
will occur at a declining rate, since they involve lessening changes in 
leaf area or canopy closure. A similar sequence of rate changes will 
occur as the plants turn from green to yellow/brown. 

The sequence of rate changes described suggests that sigmoidal 
curve shapes could be used to characterize both the greenness-increase 
and greenness-decline phases of crop spectral development. Such curve 
shapes are commonly used to describe growth and development phenomena 
in biological populations. A curve of Greenness vs. day of year for 
spring small grains is illustrated in Figure 1. 

Within the group of crops that follow this general spectral de- 
velopment pattern, significant variations in the overall temporal 
pattern can occur. The particular sequence of development stages 
through which a crop passes affects the spectral character and the 

/ 
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DAY OF YEAR 

FIGURE 1. SPRING SMALL GRAIN GREENNESS PROFILE 




temporal-spectral pattern of the crop. Variations in growing con- 
ditions can also alter the character of the individual plants and the 
canopy structure in the field, and thus again change its spectral 
qualities and temporal-spectral pattern. Some such differences can be 
adapted to by a well-chosen profile model form and can be detected by 
evaluation of model parameters. Other differences, however, may be of 
such a nature that they cannot be accurately represented by the same 
model. Different greenness measures can also behave differently in 
this regard !6], 

2.2 CP.vP CALENDAR SHIFT ESTIMATION CONCEPT 

Healthy fields of a given crop type might be expected to follow a 
common pattern of spectral development in the course of a growing sea- 
son with the result that, on any given day, spectral signals from 
nearby fields wouldn't vary to any great degree. It is more common, 
however, to encounter wide dispersion ot signal values, as depicted in 
Figure 2a. Oop calendar shift techniques are intended to account for 
that portion of signal variability which is due to differences in stage 
of development on the day of observation. As illustrated in Figure 2b, 
the highly variable data points of Figure 2a can all be fit with the 
same profile shifted al aig the time axis. By shifting all the data 
points to a common profile, as in Figure 2c, the spectral impact of 
development stage differences is adjusted for and better comparison of 
fields is made possible. 

The idea of using temporal-spectral patterns to adjust for de- 
velopment stage differences was first suggested by Dr. Cautam D. 

Ba tilt war [ 7 ] . Starting from this initial work by Badhwar, we at KK1M 
developed refinements of the crop calendar shift technique and alter- 
native approaches (8,3). First, In order to represent the actual 
pattern of spectral development of the crop of interest, we developed 
and began to use the model form (Equation 2) which is discussed in 
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Section 2.3. (Badhwar, in turn, continued development of his approach, 
adding a model form related to Equation 2 (9).) Second, we used a 
cross-correlation computation in determining the goodness-of-fi _ of a 
set of datu values to the shifted reference profile. The cross- 
correlation measure is of the* form: 


1 + 


* c , 
(£F i + , G i> 2 


(l) 


where 

Fj = reference profile function value 
s data value 
t = shift value 

This measure has the advantage of being essentially independent of 
scale and so places emphasis on the overall shape of the data profile 
rather than its amplitude which can be influenced by factors not re- 
lated to stage of development. 

Figure 3 Illustrates results obtained when the FR1M crop calendar 
shift technique was applied to actual Landsat data from spring wheat 
pixels I 10) . 


2. 3 PROF! l.K MODELING FOR SPRING SMALL GRAINS 

Any of a variety of mathematical model terms might he used to re- 
present the temporal -spectral characteristics of crops. The initial 
model form developed at ER1M for characterising green development 
profiles of spring small grains is of the following form [2,3,8,10]: 

F(t) » at* \* 1 1 



( 2 ) 
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FIGURE 3. 


RESULTS OF CROP CALENDAR SHIFT ESTIMATION APPLIED TO ACTUAL 
LANDSAT DATA (Spring Wheat Pixels) 
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where 


and 


F(t) ■ profile value at time t 

t * (day of year) - (day of first detectable 

green development) 


a,b,c - model parameters, with b and c determining 
the shape of the profile. 


This model form best exhibits the desired double-sigmoid shape for 
representing green development, as illustrated in Figure 1, when the 
spectral variable being used has values near zero at t~0. This implies, 
in general, a need for both a spectral and a temporal offset, the first 
being the value of the green measure for t < 0 (i.e., the value for 
bare soil) and the second being the day of first detectable green 
development. 


Specification of the bare soil value, i.e., the spectral offset, 
is simplified by use of data normalization to reduce variability caused 
by external (non-crop) effects (as discussed in Section 3.3) and use of 
Tasseled-Cap Greenness as the green measure because it tends to keep 
constant values in the presence of soil reflectance variations. The 
CCSHFT subroutine described in this report expects to receive, as in- 
puts, Tasseled-Cap Greenness values computed from Landsat multispectral 
scanner (MSS) band values according to the following relationship* [5]: 

G = -0.2837 MSS4 - 0.66006 MSS5 + 0.57735 MSS6 (3) 

+ 0.38833 MSS7 + 32. 


where the Landsat MSS band values preferably have been normalized 
against externally caused variations (Section 3.3). Values for the 
profile variable in CCSHFT are computed as follows: 

F a - G i - 25. (4) 


*UVCIE Landsat-2 calibration of the data is assumed. 


This spectral offset has been found satisfactory for use in analyzing 
and processing spring small grains data from the Northern U.S. Great 
Plains [2,10]. 

Determination of the temporal offset, which can vary from sample 
to sample, can be accomplished with a crop calendar shift estimation 
procedure, as previously discussed. 

The reference profile of Greenness that is built into the code of 
CCSHFT was determined through regression analysis and leaBt squares 
fitting of data values from several 5x6-mile segments to the model form 
of Equation 2. 
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CONSIDERATIONS AND LIMITATIONS 


Three general categories for consideration in the application of 
temporal-spectral characterization techniques are level of accuracy, 
model form attributes, and data normalization, 

3.1 LEVEL OF ACCURACY 

Temporal-spectral profiles may be used for a wide variety of pur- 
poses. Corresponding to this variety is a range of accuracies required 
from the model form. For some purposes it may suffice to roughly 
approximate the general shape of the profile, even to the point of 
simply connecting the observations with a sequence of straight line 
segments. At the other extreme, a highly accurate fit of physically- 
based mathematical form may be required, drawing on growth and develop- 
ment models, weather data, and other types of data in addition to the 
spectral samples. Between these extremes is a range of profile appli- 
cations with a range of accuracy demands. Each new application must 
take into account its required level of accuracy, and any model should 
be used only at the level of accuracy for which it was developed, or 
at higher levels only after careful testing. 

For crop calendar shift estimation of spring small grains, we have 
found that the day on which the maximum profile value occurs is the 
most dominant and consistent profile feature. This peak occurs just 
prior to or at the start of heading when green leaf area is at a maxi- 
mum. A good overall characterization of the profile shape is sufficient 
for estimating the peak date, and we have found the model form of 
Equation 2 to be satisfactory. 

Other crop calendar dates, such as date of planting or specific 
growth stages, might be inferred through use of this approach but only 
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insofar as their relationship to the date of peak is constant or nearly 
so. For example, a preliminary evaluation of crop calendar shift tech- 
nology applied to field reflectance measurements from soybean fields 
found a correlation between the shift estimate and row spacing; row 
spacing influenced when maximum cover occurred and, consequently, 
affected the day of maximum Greenness. Fields planted on the same 
day were shifted different amounts, and thus planting date estimates 
based on crop calendar shift would have been in error. 

3.2 MODEL FORM ATTRIBUTES 

For a model form to be of practical use in characterizing crop 
temporal-spectral profiles, the number of parameters which must be 
estimated should be kept small. In addition, the overall characteris- 
tics of the profile shape should be easily derived from the estimated 
parameters. The model must produce stable results, and be able to 
adapt to the range of shape variations likely to be encountered. Nor- 
mal variations shouldn't adversely affect the profile estimation process 
At the same time, however, and particularly in applications for which 
the accuracy requirement is high, the model must respond to those dif- 
ferences which carry the information of interest and maintain a close 
fit to the data samples. 

Some characteristics of green development profiles of spring small 
grains and Equation 2 were noted in the preceding section. During our 
prior development [2,10] of a spring wheat labeling procedure, we com- 
puted separate profiles for wheat and barley for several segments. We 
found that crop calendar shifts computed using these different profiles 
for reference usually differed by less than two days and so were able 
to use a single reference profile for shift calculations in the pro- 
cedure. For other purposes, we have computed segment-specific profiles 
based on shifts estimated with the general reference profile [2,3j. 
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3.3 DATA NORMALIZATION 


In order to insure that profile shapes fit to Landsat observations 
are influenced only by agrophysical phenomena, and to allow comparison 
of profiles over wide geographical regions, data normalization is 
highly desirable, if not essential. Corrections for sun angle, sensor 
calibration, and atmospheric haze effects, and detection of clouds, 
defective data, etc., should be carried out prior to profile charac- 
terization, particularly when profile features are to be used for crop 
Identification or condition assessment. We recommend and use ERIM’s 
spatially-varying XSTAR haze correction algorithm [10,11) for correc- 
tion of atmospheric effects and associated preprocessing algorithms 
5,10). 
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DETAILS OF CROP-CALENDAR-SHIFT SUBROUTINE CCSHFT 

This section describes our crop-calendar-shift subroutine called 
CCSHFT. User documentation is presented in Appendix A, a FORTRAN list- 
ing in Appendix B, and test cases and results in Appendix C. 

It is important to indicate here what this subroutine does, what 
it requires, what environment it will best operate in, and some things 
it does not do. 

The subroutine program does: 

(1) Operate on multidate values of the Tasseled-Cap Greenness 
variable (see requirements and preferences below), 

(2) Determine if sufficient appropriate acquisitions are 
present, and 

(3) Compute an estimated day of peak Greenness and a goodness of 
fit, assuming the observations to be from a spring small grain. 

The subroutine requires or assumes: 

(1) A minimum of three acquisitions that are separated by 17 or 
more days and occur during the growing stages of wheat — 
after emergence and before harvest. 

(2) Tasseled-Cap Greenness values with a 32-count offset having 
been added to eliminate negative values. (See Equation 3, 

Sec. 2.3) 

The preferred operating environment provides: 

Data that have been preprocessed to reduce non-crop-related 
variations due to satellite calibration, sun angle, and atmos- 
pheric haze. We use and recommend use of the ERIM spatially 
varying XSTAR algorithm and associated correction factors. 

(See Sec. 3.3) 
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The subroutine does not: 

(1) Make a detailed estimate of any crop stage other than day 
of peak Greenness (which apparently occurs just prior to 
heading at the time of peak leaf area and/or flag leaf 
development) . 

(2) Estimate the value of peak Greenness or other specific 
profile characteristics, since they are not required for 
the intended application, crop calendar shift estimation 
for spring small grains. 

The steps followed by the subroutine are as follows: 

(1) Read in data values and corresponding days of year. 

(2) Compute rough estimate of day of peak Greenness. 

(3) Accept acquisitions within a prescribed window around 
the estimated peak. 

(4) Compute a cross-correlation function between the observed 
values and corresponding values of a standard small-grains 
reference profile, as they are shifted in time past each 
other. 

(5) Choose that shift value which maximizes the cross-correlation 
function and compute the corresponding day of peak Greenness. 

(6) Output the shift value and a goodness of fit. 

The goodness- of-f it value returned by the subroutine is derived 
from the maximum cross-correlation coefficient, and typically ranges 
from 0 to 1, where 1 represents a perfect fit. However, negative 
values are possible and indicate a very poor fit. 
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APPENDIX A 

USER DOCUMENTATION OF SUBROUTINE CCSHFT 


SOUTINE: 

CCS HPT 


VERSION: 

1.0 


DATE: 

February 1980 


PROGRAHMER: 

E. Crist £.C • 


LANGUAGE: 

FORTRAN 


INTERFACE: 

A Standard FORTRAN Integer 

Function 

PURPOSE: To compute an estimate of the 

day of maximum Greenness for 


•Mil grains targets. 


CALLING SEQUENCE : 

I ROCCSHFT ( DATA , ACQDAY , NACQ , PKDAY , GF I T ) 
or 

CALL C’SilFT (DATA , ACQDAY, NACQ , PKDAY ,GFIT) 


ARGUMENTS: 



NAME 

TYPE 

DESCRIPTION 

DATA(IS) 

R 

Tasaeled-cap Greenness values 
(with 32-count offset), or 
•99.0 if flagged by SCREEN. 

ACQDAYU5) 

I 

Acquisition days (1-366) corresponding 
to values in DATA(). 

NACQ 

I 

Number of acquisitions. 

PKDAY 

I 

Output estiMte of day of peak 
Greenness, or ’O' if no peak estimate 

can be made. 

GFJT 

R 

Output measure of data fit to profile 
based on cross-correlation factor. 




RETURN CODES: 


fc- 

j 

-L_4 


CCSHFT * 0 Successful return 

1 Too few acquisitions 

2 Too few or inadequately spaced acquisitions 
(after first rough shift) 

DESCRIPTION: 

CCSHFT carries out the following sequence of steps in estimating 
the day of maximum Greenness: 

1. Check acquisition availability — there must be a minimum of 
three non-SCREENed acquisitions. If there are less than three, 
CCSHFT = 1, and return to calling program. 

2. Fit a quadratic to three acquisitions, including the one with the 
largest Greenness value. The maximum value of the quadratic is 
the first estimate of the day of peak Greenness. If the quadratic 
cannot be fit, the acquisition day corresponding to the highest 
Greenness value in DATA() is used as the first peak estimate. 

3. Line up the data to a reference Greenness profile by matching 
the peak day estimate to the peak day of the profile. The 
reference profile is illustrated in Figure A-l. 

4. Re-checlc acquisition availability — now there must be at least 
three acquisitions in the range marked 'A* in Figure 1. In 
addition, a spacing constraint is applied, such that consecutive 
acquisitions by different satellites (i.e., approximately 9-day 
spacing) are only treated as one acquisition in meeting the three 
acquisition requirement. If there are less than three adequately 
spaced acquisitions, CCSHFT ■ 2, and return to calling program. 

5. Shift the data values in time relative to the reference profile 
(±30 days) and select that shift which maximizes the cross 
correlation term. 







FIGURE A- I 



where is the reference profile value, 

Gj, is the data value, 

T is the shift value. 


6. Compute the ffnal peak day estimate by using the selected 
shift value to adjust the first peak day estimate. 






INTEGER FUNCTION CCSHFT(DATA,ACQDAV,NACO,PKDAy,GFITt 


c****< 

Co 

C* 

c* 

(,aui 

ROUT tNC TO ESTIMATE 

>•**•••******• A************* 4*4 ••*******••«* 

DAY or PEAK GREENNESS FOR GRAIN TARGETS, 

c* 




c* 

argument 

Type 

DESCRIPTION 

c* 




c* 

DAtA(15) 

R 

tasseled-cap GREENNESS values, or 

C A 



-99.0 IF FLACGtO BY SCREEN. 

c* 




c« 

acqdavus) i 

ACGUISIUON OAVS (1.366) CORRESPONDING 

c* 



TO DATA VALUES. 

C 4 




c* 

NACQ 

I 

NUMBER OF ACQUISITIONS INPUT. 

c* 




c* 

PKOAV 

I 

OUTPUT ESTIMATED DAy OF PEAK GREENNESS 

c* 



(1-366). 

c* 




c» 

SFIT 

R 

MEASURE of fit at maximum 

c* 



cross-correlation. 

c* 




c* 




c* 

NRITTEN bt 

E, P. CRIST 

c» 


ENVIRONMENTAL RESEARCH INSTITUTE OF MICHIGAN 

c» 


P.O. BOX 

8618 

c* 


ANN ARBOR 

, MICHIGAN 48107 

c* 





C* 

C* DECLARATIONS and DIMENSIONS 

c* 

INTEGER PKDAV,SDAV, SHIFT, TDIFF,TOPDAV 

c» 

INTEGER ACQDAV(15),0AV(I5),T(3> 

c* 

REAL INt,MAX,MAXR 

c* 

REAL DAtA(|S),G(3),GR(15),PRUFIL(I50) 

C* 

DATA PROF It/0. 651,0.651,0. 651, 0.651, 0.651, 

*0. 651, 0, 651,0.651,0. 651, 0.651,0,651, 0.651,0.651, 

10.651,0 651,0. 651,0.651,0.651,0.651,0.651,0.651, 

40. 65 1. 0. 651.0.65 1.0. 65 1.0. 65 1.0. 65 1.0. 65 1.0.651, 

40.651.0. 651.1.596.2.693.1.894.5.176.6.517.7.904, 

19. 323, 1ft. 760, 12. 21 6, 13. 671, 15, 1 19, 16. 553,17.965,19.348, 
t 20. 69 7, 22,000,23,265, 24.474, 25,627,26. 719, 27.746,25,706, 
429. 595, 30. 4 1 2, St, 153, 31. 8 18, 32. 405, 32. 9|5, S3. 345, 33. 696, 
133, 972, 34. 170, 34, 29?, 34, 340, 34. 315, 34. 221, 34. 058,33. 831, 
433.541, 33 , 192, 32, 787, 32, 329, 31, 821, 31, 268, 30.673,30,039, 
429.371.28.671,27.943,27.191,26.418,25,628,24.1124,24.010, 
423, 1 88, ?2. 362, 21. 533, 20. 706, 19. A82, 19,064, 18, 255, 17. 455, 
416.668,15.894, 15.136,14.394,13.671,12.967, 12.282,11,619, 




IS 

quality 






c* 


,‘i 


H0.*7T,t0.55T,*.7O0,*.lAS,#,6S5.8.l0*,7.S*T,7.U«. 
*6.652,6.2U.5.T*6,5.«00,S,025,*.670 t *,14*,A,Ol«, 
*l,T2t,l,*#t,5.l7*,2,*il,2,70J,2,*»7,2,2B7,2.l»0, 
tl,*2S,t,76«,l,6t*,I,«75,|.SA6. I, 227, 1,117,1,016, 
l0,*?J,*.«^6,0,75* t 0.6*e,0.6?2,Q.5».2,0.6fl?,0,*S6, 

SO, *11, 0,46*. 6. 552, 0,2*8,0.267, 0,21*, 0,214, 0,1*0, 

*0.170/ 

CC$HFT*A 
C* 

c* 

C* STEP tl 

c* 

C* FIT A QUADRATIC TO THE THREE HIGHEST GREENNESS VALUES IN OROER TO 
C* SET A FIRST ROUGH ESTIMATE OF THE DAV OF PEAK GRFt NNESS. 

C* 

NllSCOaO 

MAX«0 

DO 10 lA«l,NACQ 

IF (DaTA<IA),EQ,>**.) GO TO 10 
NUSfDsNUSEDM 

c* 

c* gno contains non-screened data values, offset rv as counts, 

C* OAVO CONTAINS the CORRESPONDING ACQUISITION DATS OF YEAR, 

c* 

GR(nUSED)bDAIA(IA).2S. 

OAY(NUSED)sACQDATtlA) 

c* 

C* LOCATE HIGHEST GREENNESS VALUE, 

c« 

IF (DATAdAI.lE.HAX) GO TO 10 
HAXbDaTA(IA) 

IhAXbNUSED 
10 CONTINUE 
C* 

C* THREF nOn-SCRFFNFD ACQUISITIONS ARE REQUIRED - LESS KILL 
c* NOT PRODuCF A RELIABLE ESTIMATE of THE SHIFT, IF THREE 
C* ACQUISITIONS are not available, set ccshftm AND RETURN, 
c* 

IF (NUSfD.GE.4) GO TO is 
CCSHFtm 
PKDAY sO 
GFITKfl, 

RETURN 

c* 

15 CONTINUF 

c* 

c* THO ACQUISITIONS IN ADDITION TO THE MAXIMUM ARf TO BE USED. THE 
C* CHOICE OF ACQUISITIONS DEPENDS ON THF LOCATION OF THE MAXIMUM, 

C* IDFALLT, OnE ACQUISITION ON EITHER SIDE OF THE MAXIMUM IS CHOSFN, IF 
C* THE NAXImUm IS THE FIRST OR LAST ACO'-jISITON, HOMfvER, THE TWO 
C* NEAHFST ACQUISITIONS ARE USfO, 

c* 

IF (NUSrO«1MAX)20,20,2S 


t 


SShb: 


c* 

c» 

€4 


c* 

c* 

c* 

c* 

c» 

c» 


c* 


c* 

c* 

c* 


c* 

c* 

c* 

c* 

c* 

c* 

c* 


c* 


20 IU3E*fM*x«2 
cn TU «A 

25 IP UMAx-t >30,30*35 

34 IUSF»TMAX 
60 TO 4A 

35 lUS£«!H»x.t 
40 CONTINUE 

data arrays are further subset for ease of quadratic computation. 

IIISE*TUSE-1 
DO 50 lA*l ,3 
lUSFMUSf *1 
6(;«)>GR(!USE) 

TTCA)sDAT(IUSE) 

50 CONTINUE 

COMPUTE A STRAIGHT LINE THROUGH THE FIRST AND THIRD POINTS. IF THE 
MIDDLE VALUE FALLS BFLOW THIS line. THE quadratic computation MILL 
FIND A MINIMUM INSTEAD OF A MAXIMUM. IN THIS CASE, USE the day OF 
THE MAXIMUM GREENNESS OBSERVATION AS ThE FIRST PFAK DAY ESTIMATE, 

slopes tr,(n-Gm)/(T(t)-Tm) 

INT*G(!i-T(t)*SLOPE 

CHKVAL*SLOPE*T(2)t!NT 

IF (G(2i.r.E,CHKVAL)60 TO 55 * 

TOPOAVBDAV(IMAX) 

GO TO 70 
55 CONTINUF 

COMPUTE terms of the quadratic equation. 

0IF|»t(U-T(2) 

DlF?*T(?)-TC3) 

SUMI«T(1 j*T(2I 
S"M?*t(?)»T(J) 

DIFsQisdIFIaSUMI 

0IFSQ?sDIF2*SUM2 

DIFG1*G(1)-G(2) 

DlFG2sG;2)>G(3) 

Aa(nlF2*DIFGt-DIFI*DlFG2)/(DIF2*DlFSQ|«DiFl*DIFSQ2) 

COMPUTATION OF DAY OF PEAK REQUIRES DIVISION Bv 'A' TERM. SO 
IF A»0. THF COMPUTATION CANNOT BE MADE, IN THIS CASE, AGAIN, 

THE DAY OF THE MAXIMUM GREFNNfSS OBSERVATION IS USED AS THE 
FIRST PEAK D*f ESTIMATE. 

IF (A.NF.O.) GO TO bO 
TnPnAVRDAY(IMAX) 

GO TO 70 

so Continue 



BR(0IfG?-DIFS02*AT/01F? 


2 EMM 


TOPOAvr -1 .* 8 / (!, aa)a ,5 
C* 

c* ThF MAXImUM value of the quadratic must fall In the rangf 

C* OF THE TnRfE DATA VALUfS USED IN ITS OE TERHINATION, OTHERWISE* 
C* USE THE NEANFST OF THE THRFE AS THE PEAK DAV ESTINATE. 

C» 

IF ( TnPDAY-T ( I ) )6 1 #6S,t>5 
At TOPftAV»T(II 

GO TO 70 

AS IF (TnPDAY-T(J))70. 70.66 
66 TOPnAVsrm 

70 CQNTINUF 

c* 

c* compute thf difference hftneen the rfference profile peak 
c* Day and the estimated PEAK DAY OF THF GIVEN TARGpT* 
c* 

SMIFTrUO-TOPDAY 


c* 

c* 

C* STEP 2l 
C* 

c* THF FINAI ESTIMATE of tmf day of PFAK GREENNESS IS COMPUTED 
C* RY DETERMINING THE SHIFT Al ONG THE DAY OF YEAR Axis WHICH 
C» MAXIMIZES thf f ROSS-CORREL A T ION BETwEtN THE DATA VALUES AND 
C • A REFIREnCF profile. 

C* The REFtRtNCF PROFTtF CONSISTS OF VALUFS CORRESPONDING TO 
CA A 90-DAY InTFRVAL WHICH ENCOMPASSES ESSENTIALLY THE ENTIRE 
C* GROW INC. SEASON FOR SPRING SMALL GRAINS* AND A TO-DAY * TAtL • 

c* on either sidf. for any gtvfn set of ohservatiuns. oni t those 
C* WHICH Fall In thf 90-day INTFRVAL (USING THE DAV OF PEAK 
c* GHEEnnFSS FSTIMATE JUST COMPUTED) ARE CONSIDERED IN THE 
CA CROSS-CORRFLATION CALCULATION, 

CA 

NSPRD«0 


c* 

c» 

c* 


75 

so 

as 


NUSFDiO 

DO 9Q I A> T * NACQ 

IF (DaTA(TA),EQ,-99.) GO TO 90 
SDAv>ACQDAY(IA)aSHIFT-9« 

IF fSDAV.LE.30 .OR. SDAY.GT.I20) GO TO 90 

OrSERVATIONS must RE MORE than 9 DAYS APART (SEE BELOW), 

IF (NUSED.GT.O) Gn TO 75 
Tr»IFF«IR 
Gfl TO 60 
CONTINUE 

T0IFF*SDAY-0AY( IPRf V) 

continue 

IF (TDIFF.lt, 15) GO TO 85 
IFRFVsNUSEDaI 
NSPRDsNSPRDa I 
CONTINIIF 

NMSFDaNUSEDaI 

GR(NUSED)*OATA(TA)-25. 





2 ?ri 


c* 

c* 

c* 

c* 

e« 

c* 

c* 

c» 

c* 

c* 

c* 

c* 

c* 

c* 

c* 


o*y(Nuaeo)*sn*y 

*o continue 

second CMcn< or «cauvsi Ttnw availability, requirements arei 
i) th«fc acquisitions in the ro-dav profile interval 

?) HfiMF than 9 OAV SPACING BETWEEN SUCCESSIVE ACQUISITIONS. 
THE SECOND CONSTRAINT fS the RESULT OF OBSERVATIONS THAT 
IN TERMS OF THE SHIFT CALCULATION, DATA WHICH A Hr TOO ClOSFLV 
SPACED DO mOT BEHAVE AS INDEPENDENT OBSF RVATION3, TMC CUT-OFF 
OF is OATS USED IN the PROGrAm IS NOT absolute, but has 
chosen TO FLIMINATF CONSECUTIVE ACOUlStllONS FBOm different 
SATELLITES nMIlE ALLOWING CONSECUTIVE ACOUISIIIUnS FROH THf 
SAME SATFLllTE. 

ACOUlSI T |OMS NOT HEFTING THE SPACING 
IN THE SHIFT CALCULATION, IF THERE ARE 
ACQUISITIONS. 


CRITERION ARE STILL USED 

enough heil -spaced 


100 


IF (NSRrD.GE.1) 
CCSHFt*? 
PKOAVbO 
GFITAft, 

RFTORn 

CONTINUE 


GO TO 100 


c* 

c» 

c* 

c« 

c* 


ALL SHIFTS OF 
QUADRATIC ARE 


♦ OR • 30 OATS FROH THAT 
CONSIDERED. THE MAXIMUM 


COHPUTED Hi Th THE 
CROSS-CORRELATION, 


AND THE CORRESPONDING SHIFT, ARE UPDATF0 AS APPROPRIATE, 

MAaRaO. 

ISTFPt-II 
DO 1 2* TST«|,6) 

I S TERt I STEP* I 

Sl|upH« S ^ ( 

SS.'iPW^BO, 

SS JOAts*. 

00 11- IDil.'.uSFD 
S n A*«f'A»( TtiWrSTFR 
5 AU«PRD*CR( IDT -PROF IL(SOAV) 

SS f .»ROtS5f)P®(l»P»l)r II f SnAV)APROFIL(SDAV) 

SA i-)AT«SSDOAT»GR( ID)*GRIID) 

110 C'Kfl>.i<r 


C» 

C* 


R»E, /M . aSSOD* r/SlIMPRD* ISSOPRO/SUmPRD) ) 


IF fR.LT.-AWR) GO TO l?ft 
HAWPIO 

NSHfFT«IST|P 
t?0 CONTINUE 

c* 

c» The final SHIFT IS the COMBINATION OF the FIRST f based on THE 
5 * 0 U *0»* T IC » # ANf) THE SECOND (BASED ON TmE CROSS-CORRELATION TO 
C* THF REFERCnCF PROFILE). THf f STIHATFO DAT OF P£aK GREENNESS 
C* IS THE RFFfRfNfE DAT OF PEAK GREENNESS (DAV 160) ADJUSTED RT 
C • THE FSTIMATED SHIFT, 

c» 


c* 


PKf)Ay*l fc o-(SHIFTfNSHlFT) 

GFITpIO.ahawR-9, 

rftiirn 

END 


/? 

A 7 



APPENDIX C 


TEST CASES AND TEST RESULTS FOR SUBROUTINE CCSHFT 


2p 


TEST CASES FOR ROUTINE CCSHFT 


A series of test cases was run to verify the proper functioning 
of CCSHFT. Attached are listings of the testing program 'TESTCC.F', 
test data 'TESTDATA', and output results ' TEST. RESULT' . The cases 
and their purposes are as follows: 

Case 1 - test normal flow of program 

Case 2 - test first check of acquisition availability 

(lines 74-101) 

Case 3 - test mechanism for assuring that the quadratic 

estimates a maximum rather than a minimum 
(lines 134-141) 

Case 4,5 - test selection of acquisitions for quadratic, 

and mechanism for adjusting peak estimate 
to fall in data range (lines 111-118,173-178) 

Case 6 - test second acquisition check - less than 3 acqs 

Case 7 - test second acquisition check - minimum spacing 

requirements (lines 204-244) 

Case 8,9 - illustrate GFIT values from non-typical data 
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TESTDaTa 


tftftftftft 

* * * * * * 

19 : 59 t 90 

******** 

**** 

0 ?- 2 *U 

****** 

80 
* * « 

ftftfti 

i * 

1 

1 59 

157 

175 

193 

211 





2 

95 

.0 

60,0 

55 

.0 

90 

.0 

3 ft . 

0 

3 

139 

157 

175 

193 

211 





9 

95 

. 0 - 

99,0 

55 

• A m 

99 

.0 

- 99 . 

0 

s 

139 

157 

175 

193 

?1 1 





6 

60 

.0 

95,0 

55 

.0 

9 0 

.0 

30 . 

0 

7 

139 

157 

175 

193 

211 





8 

30 

.0 

30.0 

90 

.0 

55 

.0 

65 . 

0 

9 

139 

157 

175 

193 

?11 





10 

65 

.0 

55.0 

90 

. 0 

30 

.0 

30 . 

0 

11 

139 

157 

175 

193 

211 





12 

60 

,0 • 

99,0 

90 

• 0 - 

99 

.0 

30 . 

0 

13 

139 

157 

165 

193 

?1 1 





19 

95 , 

.0 

60.0 

55 

• ft » 

99 

.0 

- 99 , 

0 

15 

139 

157 

175 

193 

211 





16 

« 5 , 

.0 

95.0 

95 , 

.0 

95 

.0 

« 5 . 

0 

17 

139 

157 

1 75 

193 

211 





18 

55 , 

.0 1 

50,0 

95 , 

,0 

50 , 

,0 

55 . 

0 






X 




Tf ST, RfSijLT 
l«:59;«ft 0?-?8-fl0 


**»tf* t *t*t****tt*t****** 


CASE 1 


RETURN 

CHDEs 

0 

PKDAYs 

161 

GFITs 

CASE 2 

RETURN 

cnoFs 

1 

PKDAYs 

0 

GFITs 

CASt 7 

RE TURN 

COOFs 

0 

PKDAYs 

15? 

GFTTs 

CASE U 

return 

COOFs 

2 

PKDAYs 

0 

GFMs 

CASE 5 

return 

COOFs 

0 

PKOAYs 

1«1 

GFTTs 

CASE 6 

return 

COOFs 

2 

PKDAYs 

ft 

GFITs 

CASE 7 

return 

COOFs 

2 

PKOAYs 

ft 

GUTs 

CASE 6 

return 

COOFs 

0 

PKDAYs 

1 60 

GFITs 


* * * 
ft • 
* * 
• * 
• ft 
ft ft 
ft ft ft 


0 .99Si|6«89 


0,0 

0.32S9S23? 


0.0 


0. 99728203 


0.0 


0.0 


ft.0«9'456?l 


CASE 9 

RETURN CODE 


0 


PKOAYs 15*5 


GFTTs 0. 29122353 


